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ABSTRACT 

in 

O 

, We consider the spatial clustering of massive black hole (MBH) mergers, and 

discuss possible ways to use gravitational wave observations in the LISA and 
DECIGO/BBO range for obtaining cosmological and cosmogonical information. 
\q . Constraints on large scale structure (LSS) and merger histories may be possible 

through the detection of an alignment of the GW polarization direction with 
! principal axes of the LSS. Constraints on the merger physics and the reionization 

^ ! epoch may be obtained by GW measurements of MBH correlation lengths, in 

CO '. the case when the MBH angular momentum loss occurs through gas drag. Such 



O 

in 
c3 



c3 



measurements would provide information about the LSS and the reionization 



epoch, as well as about the astrophysics of MBH mergers, additional to and 
q . independent of that obtained from electromagnetic signals. 

Or 



Subject headings: black hole physics — cosmology: theory — galaxies: evolution 
— galaxies: nuclei — gravitation — quasars: general 



^ | 1. Introduction and summary 



Gravitational waves (GWs) will provide new eyes for studying the universe in the 21st 
century. Ground-based laser interferometers such as TAMA300 and LIGO have already be- 
gun operations in the 10 Hz - kHz band, while the Laser Interferometer Space Antenna 
(LISA) 3 , covering the 10~ 4 - 10~ 2 Hz band, will be launched in 2011. Other recently pro- 
posed facilities are the Decihertz Interferometer Gravitational Wave Observatory (DECIGO) 
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(Seto, Kawamura, & Nakamura 2001) and the Big Bang Observer (BBO) 4 , operating around 
decihertz frequencies (10 -2 - 10 Hz). These detectors may determine the spatial distribu- 
tion of GW sources, which should contain important cosmological as well as astrophysical 
information (e.g., Ioka, Tanaka, & Nakamura 2000; Seto & Cooray 2005). 

In this paper we propose novel ways to use a 3D map of GW sources, in particular 
of Massive Black Hole (MBH) mergers. MBH mergers are one of the most powerful GW 
sources in the universe. LISA may detect MBH mergers involving masses 10 3 ~ 10 7 M Q 
out to high redshifts, while DECIGO/BBO may localize 1 ~ 1O 3 M MBH mergers with an 
arcmin resolution (Takahashi & Nakamura 2003). They can also measure distances to MBH 
binaries through the change of the GW frequency during the observation. 

The MBH merger history has large uncertainties derived from the seed MBH formation 
and the angular momentum evolution of MBH binaries. The presence of high redshift quasars 
may imply that seed MBHs were formed much before the cosmological reionization epoch 
as Population III (Pop. Ill) remnants (Madau & Rees 2001; Volonteri, Haardt, & Madau 
2003). In such a scenario the event rate could reach ~ 10 4 events per yr (§ 2; Islam, Taylor, 
& Silk 2004), which makes them potentially interesting for doing cosmology. Although we 
will find it impossible to measure the spatial correlation of host galaxies of MBH mergers 
through GWs alone (§ 2), we may still use the 3D maps of MBH mergers to study the MBH 
evolution. 

One possible way to extract information from a 3D map of GW sources is the cross- 
correlation of GWs with electromagnetic (EM) sources (§ 3). We will suggest that the orbital 
axis of the MBH binary can be aligned with a principal axis of the large scale structure (LSS) 
in which the binary resides, and that this may be detectable with LISA and DECIGO/BBO. 
Such a detection would provide useful clues to the angular momentum gains of MBH binaries. 

GW sources could have a larger correlation length than that of host galaxies, and in 
this case we may detect a spatial correlation of GW sources (§ 4). We will suggest that the 
MBH merger rate may drop inside an ionized bubble, since ionizing photons heat the gas 
which could be responsible for the angular momentum loss of MBH binaries. Since ionized 
bubbles may be larger than the LSS scales at the end of reionization, we could detect the 
spatial correlation of MBH mergers and thus probe the effectiveness of MBH binary angular 
momentum loss via gas. The merger rate may also drop after reionization, thus providing 
an independent marker for this important epoch. 
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2. Event rate and correlation length 



We first estimate the event rate of MBH mergers, adopting throughout a ACDM cos- 
mology with {Q m , n A , Q b , h, a 8 ) = (0.27, 0.73, 0.044, 0.71, 0.84) (Spergel et al. 2003). Since 
little is known about the MBH merger history, we concentrate on a model in which MBH 
mergers are associated with the mergers of host halos of dark matter (Begelman, Blandford, 
& Rees 1980; Volonteri, Haardt, & Madau 2003). By using the extended Press & Schechter 
(1974) (PS) formalism, we can estimate the halo merger rate as in equation (4) (Lacey & 
Cole 1994; Haehnelt 1994). 



The comoving number density of halos in the mass range M 
given by the PS mass function as 



M + dM at redshift z is 
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where p m = 3iJof2 m /87rG ~ 2.78 x lO n f2 m /i 2 M Mpc~ 3 is the mean comoving matter 
density, a = a(M) is the rms mass fluctuation on a mass scale M at the present epoch, 
5 c (z) = 5 c D(z = 0)/D(z) is the critical linear overdensity evaluated at present for a spherical 
perturbation to collapse at z, and D(z) is the linear growth rate. These halos of mass M 
are incorporated into a larger halo of mass M 2 ~ M 2 + dM 2 at redshift z at a rate given by 
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where a 2 = a(M 2 ) (Lacey & Cole 1994). Therefore the merger rate of halos of mass larger 
than M m i n with halos of the same mass or larger per comoving volume is given by 
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This yields the total number of merger events per observed time t Q b s per redshift z as 



dN, 



merg 



dzdt, 



obs 



1 dn mcrg dV c 
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(3) 



(4) 



where dV c /dz = Ancd 2 L \dt / dz\/ (l+z) is the comoving volume per redshift, d^ = c(l+z) / Z (1+ 
z)\dt/dz\dz is the luminosity distance, and \dtfdz\~ 1 = (1 + z)H y/Q m (l + z) 3 + Qa- We 
follow Kitayama & Suto (1996) to estimate 5 C and a(M). 

We note that the merger rate in equation (3) is dominated by equal mass halos (M 2 ~ 
2M). Thus the rate is not reduced much even if dynamical friction is inefficient in initially 
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bringing halos of small mass ratios together (Lacey & Cole 1994; Yu 2002). The dynamical 
friction also becomes ineffective as the MBH binary hardens. It is unclear how the MBH 
binary loses its angular momentum afterward. The two main possibilities are through inter- 
actions with stars or with gas. Stars may not be effective since the MBH binary will quickly 
eject all stars on intersecting orbits (Milosavljevic & Merritt 2003; Merritt & Milosavljevic 
2004), while gas may reduce the binary separation efficiently (Gould & Rix 2000; Armitage 
& Natarajan 2002; Escala et al. 2004). Here we assume that two MBHs coalesce soon after 
the host halos merge. In § 4 we will propose observational methods to test whether gas is 
responsible for the angular momentum losses. 

The event rate in equation (3) is dominated by halos with the minimum mass M ~ M min . 
It is not well known what fraction of galaxies (and down to what mass) harbor MBHs, 
although almost all nearby massive galaxies have central MBHs (Magorrian et al. 1998). 
Here we assume that all halos of mass above M min (z) have a central MBH and consider three 
possibilities for M min (z): (a) The minimum halos have a virial temperature T vir ~ 10 4 K, 
above which baryons can cool via atomic hydrogen lines (Barkana & Loeb 2001; Menou, 
Haiman, & Narayanan 2001). We calculate T v ; r according to Kitayama & Suto (1996). 
(b) The minimum halos have T vir ~ 10 3 K, above which molecular hydrogen (H 2 ) cooling is 
possible (Barkana & Loeb 2001; Bromm & Larson 2004). (c) The minimum halos contain 
more than one Pop. Ill MBH on average. We assume that one MBH forms in each 3a halo 
at z = 20, and the mean number of MBHs in a given halo is calculated by equation (3) of 
Madau & Rees (2001). 

Figure 1 shows the number of merger events per yr per unit z as a function of z. We 
can see that in total ~ 10 3 events/yr are expected in case (a), and ~ 10 4 events/yr in the 
case (c). In the case (b) we plot only the region of z > 15 with a solid line because H 2 
is fragile and can be destroyed by photons in the Lyman- Werner (11.2-13.6 eV) bands well 
before the cosmological reionization (Haiman, Abel, & Rees 2000; Omukai & Nishi 1999). 
The H 2 history is very uncertain because free electrons can promote H2 formation (Bromm 
& Larson 2004), so that > 10 3 events/yr at z > 15 is not implausible. The mass density 
of MBHs in the case (b) exceeds that found in the nuclei of nearby galaxies (Madau & Rees 
2001; Yu & Tremine 2002). Thus in the case (b) the merging efficiency should drop at low 
redshift and many MBHs are to be wandering without being observed. The event rate in 
the case (c) is larger than that in the case (a). This suggests that the MBH merger might 
happen in nonluminous dwarf galaxies, which have not been explored yet. 

The MBH mass, especially in small halos, is quite uncertain. If we extrapolate the 
relation of MBH and halo mass (Ferrarese 2002), we have M MBH ~ 10 7 M Q (T vir /10 6 K) a 
where a ranges from ~ 5/2 (Merritt & Ferrarese 2001) to ~ 2 (Gebhardt et al. 2000). For 
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a ~ 2 the MBH mass in halos of T vir ~ 10 4 K is ~ 1O 3 M , so that almost all events of case (a) 
may be detectable by LISA (e.g., Hughes 2002). For a ~ 5/2, we have M MB h ~ 1O 2 M . This 
is below the sensitivity of LISA but above that of DECIGO/BBO (Takahashi & Nakamura 
2003). In the case (c), the Pop. Ill MBH mass is also unclear. Depending on the gas 
accretion onto a protostellar core, the final mass of a Pop. Ill star and its end product may 
reach ~ 1O 3 M (Omukai & Palla 2003). In this case LISA may detect almost all Pop. Ill 
MBH mergers. 

The first question to be addressed is whether the correlation length of the MBH mergers 
is measurable or not. Since the merger rate is dominated by the minimum halos of mass 
M min , the correlation length r of the MBH mergers is roughly that of the minimum halos. 
We can then estimate the correlation length r by 5 

6(M min , z)a{r )D(z)/D(z = 0) = 1, (5) 

where b{M, z) = 1 + 5 C _1 [5 2 (z)a~ 2 — 1] is the bias parameter for halos of mass M at redshift 
z (Mo & White 1996). Figure 1 shows the correlation length r as a function of z. Whether 
the correlation length r is measurable or not is determined by the number N v& i x of merger 
pairs whose distances are less than r$. We need at least N v ^ r > 10 for a 3<r detection since 

— 1/2 

the Poisson error of the correlation function is given by A£/£ ~ (Peebles 1980, § 48). 

Unfortunately we can show that Ap air at each redshift bin Sz ~ 1 is less than unity even if 
we observe for 10 yrs (even for the case (b)). Therefore we conclude that it is impossible to 
measure the correlation length ro of the host halos, using the MBH merger GWs, regardless 
of the detector sensitivity. 

However, there are other types of spatial information and other means to extract them 
involving MBH merger GWs. One of these is the exploitation of cross-correlations between 
GWs with EM sources (§ 3). Another is in the case when the spatial distribution of MBH 
mergers has a larger characteristic scale than the correlation length r of the associated halos 
(§ 4, see also Khlopov, Rubin, & Sakharov (2005)). We discuss these two cases below. 



3. Orbital alignment with large scale structures 

Among various ways to cross-correlate GWs with EM sources, the one involving the 
GW polarization and the EM LSS may be interesting. In the hierarchical structure forma- 
tion scenario, sheet- or filament-like LSSs are initially formed. Such LSSs have a preferred 
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direction, which could imprint a preferred orientation of the MBH binary relative to the LSS 
(see Figure 2). The orientation of the MBH binary may be determined through the GW 
polarization (Cutler 1998), while the associated LSS can be mapped with wide-field cameras 
such as Suprime-Cam on the Subaru Telescope (Kodama et al. 2001; Ouchi et al. 2005). 
This offers the prospect of correlating the MBH binary axis with the direction of the LSS. 

The number of events needed to detect a correlation may be small. Let Q{ (0° < Q{ < 
90°, % — 1, • • • , N) be an angle between the binary axis and the associated filament direction. 
(We focus on filaments to be concrete.) The anisotropy may then be quantified by (Struble 
& Peebles 1985) 



For an isotropic distribution, 5 ~ and the standard deviation would be ag = 90/(12iV) 1//2 
deg. Thus only N ~ 10 events may be enough for a 5er detection if the correlation is strong 
(9i ~ 0° or 90°). However a 100% alignment would be unlikely. If the alignment is ~ 20% 
(which corresponds 5 ~ 10) as suggested for the alignment between galaxies and the LSS 
(Pimbblet 2005; West, Jones, & Forman 1995; Lee & Pen 2002), we need N ~ 200 events 
for a 5o" detection. Note that the alignment between the radio jets and the galaxy disk has 
not been measured, but this is still consistent with the ~ 20% alignment because of the 
Poisson error (Kinney et al. 2000). 

We may apply the discussions on the alignment of galaxy spin axes with LSSs (e.g., 
Lee & Pen 2002; Navarro, Abadi, & Steinmetz 2004; West 1994). In hierarchical models 
of galaxy formation, the galactic angular momentum is produced by tidal torques from the 
surrounding matter in the early protogalactic stage (Peebles 1969; Doroshkevich 1970; 
White 1984). The tidal torque theory predicts that galaxy spins tend to align with the 
intermediate principal axis of the LSS, suggesting that the MBH binary may also favor an 
orbital axis normal to the LSS filaments 9i ~ 90°. Alternatively, the binary may have Q, L ~ 0° 
if galaxies merge in the direction of the collapse from a sheet to a filament. Therefore the 
detection of such correlations would be useful for probing the MBH merger history. 

A typical comoving scale of LSS at z ~ 1 is ~ 20 Mpc, which is about ~ 1° on the sky 
and Az ~ 0.01 at z ~ 1. Although LISA may not attain such precision in the localization 
of MBH mergers of mass ~ 1O 3 M (Hughes 2002; Cutler 1998; Seto 2004), we may be 
able to identify associated LSSs if we can find EM counterparts to the MBH merger (Tipler 
1975; Milosavljevic & Phinney 2003). On the other hand DECIGO/BBO should be able to 
identify associated LSSs because of its arcmin resolution (Takahashi & Nakamura 2003). We 
may determine the orientation of the LSSs at the position of the MBH merger with wide-field 
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telescopes one by one, or we may have already known the LSSs from future high-redshift 
surveys or 21 cm surveys. 

4. Effect of photoionization on MBH-drag 

One of the major problems in the MBH evolution is uncertainty about how the MBH 
binary loses its angular momentum. The two main candidates for extracting the angular 
momentum are stars and gas. With stars there is the so called final parsec problem, in 
which the binary quickly ejects all stars on intersecting orbits and cuts off the supply of 
stars (Milosavljevic & Merritt 2003; Merritt & Milosavljevic 2004), although this regions 
may be refilled via star-star encounters (Milosavljevic & Merritt 2003) or chaotic orbits 
(Merritt & Poon 2004). On the other hand, gas may be efficient in reducing the binary 
separation (Gould & Rix 2000; Armitage & Natarajan 2002; Escala et al. 2004). Here 
we assume that gas is responsible for the angular momentum losses, and show that this 
assumption may be verified or disproved by the spatial distribution of MBH mergers. 

We point out that a photoionizing background would heat the gas that is assumed to be 
responsible for the angular momentum losses and would inhibit the drag exerted on MBHs 
by gas, especially in small halos. A photoionizing background photoevaporates gas in mini 
halos with virial temperatures T vir < 10 4 K and substantially reduces the cool gas in halos 
with 10 4 K < T vir < 10 5 K (Ikeuchi 1986; Efstathiou 1992; Thoul & Weinberg 1996). 
Since the amount of MBH-dragging gas is reduced, the MBH merger rate should drop in an 
ionized bubble, possibly to almost zero because the event rate is dominated by small halos. 
This would result in a spatial correlation of MBH mergers, and its detection would indicate 
that MBH binaries lose their angular momentum via gas. 

The evolution of ionized bubbles during cosmological reionization is still unclear, despite 
intensive studies (Barkana & Loeb 2001). The analysis of Lya spectra in the highest redshift 
quasars indicates that the reionization ends at z ~ 6 (Fan et al. 2002), while the WMAP 
polarization data imply an onset of reionization at a much higher redshift z ~ 17 ±5 (Kogut 
et al. 2003). Also under discussion is whether the reionization proceeds from high- to low- 
density regions (Furlanetto, Zaldarriaga, & Hernquist 2004) or from low- to high-density 
regions (Miralda-Escude, Haehnelt, & Rees 2000). 

Recent observations and simulations suggest that the ionized bubble size at the end 
of their overlap is essentially determined by the finite bubble size and the cosmic variance 
(Wyithe & Loeb 2004; Furlanetto, Zaldarriaga, & Hernquist 2004). A comoving size R of 
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a finite bubble occupies a redshift interval of 

dz 



R 

Az 



c(l + z) dt ' (7) 

which produces a scatter ~ Az of reionization redshift. On the other hand, reionization in 
a certain region would be completed when the fraction of collapsed mass exceeds a critical 
value in this region. Because of the cosmic variance, a different region reaches the critical 
point at different redshift with a scatter 



.2 



. -l + A /l--iL (8) 
1 + 2 <J c (z) V <jf 



where or = <r[M(i?)] is the rms mass fluctuation over the comoving radius R at the present 
epoch, the mass M(Rq) within Rq is the minimum galaxy mass and Sr is the mean overdensity 
on the radius R. Here we set 5r = a R , and choose M(R ) to have T vir = 10 4 K. 

Comparing equation (7) with equation (8), we find that the scatter in the reionization 
redshift is Az ~ 0.14, and the ionized bubble has a comoving size R ~ 60 Mpc (9 b ~ 0.4° on 
the sky) at the end of reionization at z ~ 6. This scale is much larger than the correlation 
length of host halos r , and may be detectable through GWs. Let us consider here a redshift 
shell z ~ z + Az at the end of reionization (see Figure 3). Since few mergers take place inside 
an ionized bubble and the bubble separation is about their size R ~ 60 Mpc (6 b ~ 0.4° on 
the sky), the MBH mergers should have a correlation £ ~ 1 on this scale. The number of 
merger pairs whose distances are less than R (less than 9j, on the sky) is approximately 

2 f _ 20 f dN mcrg /dzdt ohs y /T obs y / Az_y 2 

4 ~ V W 3 /yr/z J V 10 y r / V - 14 / V°- 4c 



^Vpair ~ 



dN r 
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dzdt, 



mCrg T A y 
1 nhelAZ 



obs 



&t_ ^ 20 ( dN mcrg /dzdt ohs \ 2 / T obs \ 2 / Az \ 2 / jh_ i , , Vi 

4 ~ V w 3 /yr/z J V 10 y r / V - 14 / V - 4 °/ 

for an observation time T b s - Since the Poisson error of the correlation function is A£/£ ~ 

— 1/2 

A/pair , the spatial correlation of MBH mergers may be detectable by GWs at a 4<r level in the 
case (b) of Figure 1. (Even if the actual bubble size is less than ~ 60 Mpc (Miralda-Escude, 
Haehnelt, & Rees 2000; Gnedin 2000), we may cross-correlate the merger position with the 
bubble position to detect a correlation as discussed in § 5.) To locate MBH mergers within an 
angle 6 b LISA may need EM counterparts to the MBH mergers, while DECIGO/BBO, with 
a resolution better than 9 b (Takahashi & Nakamura 2003), may do so without counterparts. 

The bubble overlap may also be accompanied by a sharp drop in the merger rate because 
the photoionization effectively raises the minimum halo mass that harbors a central MBH 
(e.g., from the case (b) to the case (a) in Figure 1). Therefore the merger history may also 
be used to determine the reionization redshift. 6 



6 Wyithe & Loeb (2003) also suggest that the merger rate drops after reionization. However their reason 
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5. Discussion 

A precise estimate of the reduction of the merger rate at reionization is beyond the scope 
of the paper, since it would depend on the MBH-dragging mechanism (Escala et al. 2004) 
and the degree of self-shielding of ionizing photons (Susa & Umemura 2004). We note that 
gas will eventually evaporate even if the gas is self-shielding in mini halos with T vir < 10 4 K. 
This is because the outside layer of gas is always ionized under the ionizing background and 
can escape from halos. We also note that the photoionization effect may be weak at high 
redshift z > 10 (Dijkstra et al. 2004). 

In the future, the spatial distribution of ionized bubbles will be determined by other 
methods such as 21 cm tomography (Tozzi et al. 2000; Ioka & Meszaros 2005), Lya spec- 
troscopy (Miralda-Escude 1998) and dispersion measures (Ioka 2003). For such measure- 
ments, the cross-correlation of the merger position with the bubble position (like the cross- 
correlation with LSSs outlined in § 3) would make it easy to detect a signal. For example, if 
the fraction of ionized bubbles on the sky is /& ~ 1/3 in a redshift shell z ~ z + Az, we only 
need ~ 5 2 / b _1 ~ 80 events in z ~ z + Az to detect a 5a correlation. The rate of MBH merg- 
ers and the mechanism by which the angular momentum loss occurs may be independently 
constrained or determined through purely GW measurements, such as outlined in §4 

In summary, we have discussed two possible ways to exploit gravitational wave obser- 
vations in the 10~ 4 — 10 Hz range for gleaning information about the high redshift LSS 
formation and about the mechanism of angular momentum loss in MBH mergers, as well as 
the reionization of the universe. The former would be indicated by an alignment of the GW 
polarization direction with a principal axis of the LSS, due to a preferential vector along 
which the galaxies approach. The latter would be detectable through an MBH correlation 
length larger than that of the host LSS, which can occur if the MBH angular momentum loss 
occur through gas drag. In both cases we would gain information about the LSS additional to 
and independent of any obtained from electromagnetic signals, and in the latter also about 
the astrophysics of MBH mergers as well as an independent measure of the reionization 
epoch. 

We thank M. J. Rees for useful comments. This work was supported in part by the 
Eberly Research Funds of Penn State and by the Center for Gravitational Wave Physics 
under grant PHY-01-14375 (KI), and NASA NAG5-13286, NSF AST 0307376 (PM). 



is not the suppression of angular momentum losses but the suppression of MBH formation. 
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Fig. 1. — Upper panel - Number of MBH merger events per yr per unit redshift z as a 
function of z. Three cases are considered for the minimum halos that harbor a central MBH: 
(a) The virial temperature T vir of minimum halos is 10 4 K, above which atomic cooling 
is effective, (b) The minimum halos have T vir = 10 3 K, above which molecular cooling is 
possible. This is shown with a dashed line for z < 15, since molecular cooling may be 
ineffective in this range, (c) The minimum halos contain more than one Population III MBH 
on average, where we assume that one MBH forms in each 3a halo at z — 20. Middle panel 
- Correlation length of the MBH mergers as a function of z. Case (b) is plotted only for 
z > 14, and case (c) only for z < 18 (see text). Lower panel - Minimum halo mass that 
harbors a central MBH M min (z) as a function of z. 
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Fig. 2. — Schematic figure of LSSs (hatched region) and MBH binaries, illustrating the 
preference for MBH binary axes to be normal to the filaments/sheets of LSSs. The sizes of 
MBH binaries are magnified to make them clear. 
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Fig. 3. — Upper panel - Schematic figure of reionized bubbles, with neutral regions shown 
as shaded zones. Lower panel - Schematic sky view at the end of the bubble overlap era. 
Few MBH mergers take place inside an ionized bubble if the MBH binaries lose their angular 
momentum via gas drag, as the photoionizing background reduces the amount of MBH- 
dragging gas. 



